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ABSTRACT 

We describe a wavelet transform analysis of the ROSAT PSPC images of the Coma cluster. On small 
scales, < 1', the wavelet analysis shows substructure dominated by two extended sources surrounding 
the two brightest cluster galaxies NGC 4874 and NGC 4889. On slightly larger scales, ~ 2', the wavelet 
analysis reveals a filament of X-ray emission originating near the cluster center, curving to the south 
and east for ~ 25' in the direction of the galaxy NGC 4911, and ending near the galaxy NGC 4921. 
These results extend earlier ROSAT observations and further indicate the complex nature of the cluster 
core. We consider two possible explanations for the production of the filamentary feature as arising from 
interactions of the main cluster with a merging group. The feature could arise from either ram pressure 
stripped gas or a dark matter perturbation of tidally stripped material. 

Subject headings: galaxies: clusters — galaxies: ISM — X-rays: galaxies 
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1. INTRODUCTION 

Clusters of galaxies have relatively long dynamical time- 
scales, comparable to the age of universe, and hence can 
provide insights into the conditions and properties of the 
universe from which these massive systems arose. The 
study of substructure in clusters has been applied to de- 
riving cosmological parameters (Richstone, Loeb & Turner 
1992; Kaufman & White 1993; Lacey & Cole 1993; Mohr 
et al. 1995). In principle, statistical studies of substructure 
can provide strong constraints on the density of the Uni- 
verse, although uncertainties remain including the time- 
scale for substructure to be erased and the nature of 
the initial mass fluctuation spectrum (Kaufman & White 
1993). 

The Coma cluster, as the second X-ray brightest clus- 
ter and one of the nearest rich clusters, can provide in- 
sights into the reliability of estimates of fl from substruc- 
ture studies. Long-considered as the typical dynamically 
relaxed system, X-ray studies have shown its remarkable 
complexity (Briel et al. 1992; White, Briel & Henry 1993; 
Vikhlinin, Forman & Jones 1994). Optical indications of 
substructure come from the work of Fitchett & Webster 
(1987) who demonstrated the presence of substructure and 
Mellier et al. (1988) who were able to derive the mass to 
light ratios for the core substructures surrounding the two 
dominant core galaxies. Most recently, Colless and Dunn 
(1995) used a sample of 552 galaxy redshifts to further 
clarify the merger history of the Coma cluster. Biviano 
et al. (1996) combined ROSAT and extensive optical red- 
shift and photometric samples to identify the main body 
of the cluster and suggested that the main cluster body 
was undergoing rotation. They also concluded, based on 
the difference in velocity between the NGC 4889 and NGC 
4874 groups and the cluster mean, that these groups only 
recently arrived within the cluster core. 

We have revisited the extensive ROSAT PSPC archival 
data for the Coma cluster using a new technique, particu- 

1 The so-called "long term enhancement" (see explanations in Snowden 
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larly well-suited to detecting small scale structures embed- 
ded within larger scale features. This technique, "wavelet 
transform decomposition" , allows simultaneous study of 
the image at all spatial (angular) scales and allows identifi- 
cation of the scale appropriate to the detected features. In 
this Letter, we study the small scale structure of the Coma 
image with this new wavelet transform technique. We de- 
tect a long, narrow X-ray feature and discuss possible ex- 
planations as either gas which was ram pressure stripped 
from a small galaxy group during a recent merger/collision 
with the cluster or a dark matter perturbation produced 
by tidally stripped material during a group merger. 

2. OBSERVATIONS AND ANALYSIS 

We analyzed two ROSAT PSPC pointings of the Coma 
cluster core (the overall X-ray morphology of the Coma 
cluster seen from pointed ROSAT observations is discussed 
in White et al. 1993). Images were reduced using the pre- 
scriptions and code of Snowden et al. (1994). This in- 
cludes the elimination of high background intervals, which 
are primarily due to intense scattered solar X-rays, sub- 
traction of the particle and other non-cosmic background 
components 1 , and creation of the exposure maps appro- 
priate for the given energy band which are used for flat- 
fielding and merging the images obtained from different 
pointings. The cleaning resulted in 18 and 13 ksec expo- 
sures in individual pointings. The merged PSPC image 
(top- left panel in Fig 1) shows an azimuthally asymmetric 
cluster (Briel at al. 1992, White et al. 1993) whose core 
structure we investigate in detail in this Letter. 

2.1. Wavelet Decomposition 

The technique we applied to the ROSAT PSPC 
"cleaned" images employs a set of a trous wavelet trans- 
forms, closely following the approach of Starck et al. 
(1995 and references therein) with some important changes 

et al.) and the residual solar scattered X-rays 
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in the iteration scheme. An application of the multi- 
resolution wavelet-based image filtering for the X-ray im- 
ages also can be found in Slezak et al. (1994). 

The a trous wavelet transform of scale a can be (very 
roughly) characterized as the convolution of an image with 
the function equal to the difference between two Gaus- 
sians, the first positive one with a = a, and the second 
negative one with a — 2a. The basic idea of the wavelet- 
based image filtering is to perform the set of wavelet trans- 
forms with scales a, 2a, 4a..., retain only significant wavelet 
coefficients at each scale (thus filtering out the noise), and 
perform the image reconstruction from the wavelet planes. 
For the a trous wavelet transform, the reconstruction is 
performed by simple summation of all wavelet planes. At 
each scale we retain the wavelet coefficients which exceed 
some critical amplitude which corresponds to the given 
number of standard deviations (3a in our application). 
Since the wavelet transform is linear, the noise level in 
each wavelet plane is proportional to the noise level in 
the original image. Therefore, we can derive the critical 
amplitudes for each scale using the following approach. 
We first apply the wavelet transform to a simulated image 
containing Gaussian noise with dispersion 1 and determine 
the rms variations in the wavelet-convolved images. We 
then multiply these values by the local noise in the actual 
image, which is estimated using the local background and 
exposure. Thus, the faintest detectable structures have 
uniform amplitude in terms of significance, as a function 
of both position and scale. 

To assure the most complete separation of image fea- 
tures in terms of their scale, we apply the following it- 
erative procedure. We start from the smallest scale, 
a = 30", which roughly corresponds to the FWHM of the 
ROSAT PSPC PSF. We calculate the wavelet transform 
on this first scale only, find significant features, and sub- 
tract them from the original image (recall that the back- 
transformation for the a trous wavelet transform is simply 
the summation of the wavelet planes). The next iteration 
is applied to the residual and this process is repeated until 
no new significant features are found at this scale. We then 
start the process with the scale 2a and so on. The sub- 
traction of the features detected at smaller scales before 
applying the wavelet transform at larger scales results in a 
sufficient enhancement of the basic property of the wavelet 
transform analysis — its ability to separate image features 
in terms of their characteristic scale. In particular, we are 
able to isolate even extremely bright point sources primar- 
ily at the two smaller scales (i.e. < 60"), whereas without 
subtraction, they appear at all scales and often saturate 
nearby and relatively faint extended structures. 

The process described above provides a separation of 
the original image into a set of images containing only fea- 
tures of scale a, 2a, ... significant at least at the 3cr level, 
which can be either summed together to produce an adap- 
tively smoothed image or viewed individually (or in any 
desired combination) to examine the significant structures 
at desired scales. 

2.2. Small Scale Structure of the Coma Cluster 



Contour maps of the images generated for small scales 
are shown in the bottom panels of Fig 1 superposed on an 
optical image. On the smallest scales, 60" (bottom-left), 
the contours show the extended sources surrounding the 
two brightest cluster galaxies NGC 4874 and NGC 4889 
first detected in X-rays by White et al. (1993; see also 
Vikhlinin et al. 1994). On larger scales, 120" (bottom- 
right), the wavelet analysis reveals a filamentary X-ray 
emission feature originating near the cluster center, curv- 
ing to the south and east for ~ 25' towards the galaxy 
NGC 4911, and ending near the galaxy NGC 4921. The 
surface brightness enhancement near NGC 4911 was first 
noted by White et al. (1993). 

Within the outer contour of the filament, the surface 
brightness ranges from 5-10% of the underlying diffuse 
emission from Coma itself. The formal significance along 
the axis of the filament varies from 4.5cr to 6o\ Features 
with this significance never appear in simulated Coma 
cluster images. The counts associated with the emission 
filament (within the outer contour shown in bottom-right 
panel of Fig 1) correspond to a count rate of 0.12 cts s -1 
in the energy band 0.5-2 keV. The corresponding luminos- 
ity, for Hq = 50 km sec -1 Mpc -1 and a sufficiently hard 
spectrum (see below) is 3 x 10 42 ergs sec -1 . At the Coma 
distance (z = 0.0235), the tail is approximately 1 Mpc 
(25') in length with a diameter of 160 kpc (4'). 

It is interesting to constrain the spectral properties of 
the filamentary structure. It is faint compared to the un- 
derlying cluster emission, which therefore cannot be ne- 
glected and must be carefully subtracted. We use the fol- 
lowing technique for the extraction of the filament's spec- 
trum. We perform the wavelet decomposition of the broad 
band image as described above and record all the transfor- 
mations which have been applied. The identical transfor- 
mations are performed to the images in individual energy 
bands defined by the energies 0.2-0.4-0.5-0.7-0.9-1.3-2 
keV. Since the wavelet transform is linear, the relative 
fluxes in different energy bands are conserved and can be 
used to fit the spectrum. Unfortunately, there is a seri- 
ous problem with the PSPC spectral calibration for these 
particular observations, which results in anomalously low 
temperatures for the total cluster emission. 2 We therefore 
use the ratio of the filament brightness to the underlying 
cluster brightness as a function of energy to constrain the 
filament spectrum. This ratio is constant within the errors, 
which suggests that the filament has the same spectrum 
as the Coma cluster. Quantitatively, a 95% lower limit 
to the temperature in the filament is 4 keV, assuming a 
cluster core temperature of 8 keV (Hughes, Gorenstein & 
Fabricant 1988). 

3. DISCUSSION 

The wavelet analysis of the ROSAT PSPC X-ray im- 
ages of the Coma cluster demonstrates the existence of a 
new feature possibly associated with dynamical activity 
in the cluster core. We can estimate the physical parame- 
ters associated with this filament of emission assuming it 
arises from hot gas. Let us assume that we are viewing the 
filament in the plane of the sky and it is a uniform cylin- 

2 To our knowledge, a similar effect is found for other clusters observed around July 1991: A3558 (Markcvitch & Vikhlinin 1997) and A1795 
(L. David, private communication), have anomalously low best fit temperatures. This suggests that a likely explanation is a detector gain 
variation around July 1991, which is not accounted for by the present versions of the PSPC calibration. 
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Fig. 1. — (top-left) The X-ray image of the Coma cluster in the energy band 0.5-2 keV, smoothed with the Gaussian a = 1', over- 
layed on the optical image, (bottom-left) Structure detected in the Coma cluster core on the wavelet scale 60" The most prominent 
features are extended sources around NGC 4889 and NGC 4874, as well as a few point sources, (bottom-right) Structure detected 
on the 120" scale. The tail of enhanced emission originates near NGC 4874, goes to the South and East towards NGC 4911 and 
ends near NGC 4921. The formal statistical significance along the axis of the tail varies from 4.5a to 6a, and the surface brightness 
is ~ 10% of the underlying cluster emission, (top-right) Small scale structures shown in bottom panels were subtracted from the 
Coma cluster image, and the result was smoothed with the Gaussian a = 1'. After the subtraction of the small-scale structure the 
cluster core appears much more regular. The complex morphology (elongation, and the rotation of the elongation axis) therefore 
reflects the "hidden" small-scale structure, not the state of the main cluster body. 



der of radius 100 kpc and length 1 Mpc. The observed 
10% surface brightness enhancement implies that the den- 
sity of the filament must exceed the cluster gas density at 
least by about 5%, if the gas is not particularly cool (T > 
a few keV). The relative brightness of the filament does 
not change much along its axis, therefore the gas density 
enhancement in the filament is an almost constant small 
perturbation and one can use the main cluster density pro- 
file (Briel et al. 1992) to estimate the filament's gas mass. 
We derive a total average gas density in the volume de- 
fined by the filament of approximately 10 -3 cm~ 3 with a 
range (along the filament's axis) from about 2.5 x 10~ 4 to 
2.5 x 10~ 3 cm~ 3 . This density translates into a gas mass 



of 5 x 10 11 Mq. We consider two possible scenarios for the 
origin of this filamentary feature. 

First, let us suppose that this feature arises from gas 
that has been ram pressure stripped from a merging group. 
The observed gas mass is well within the reasonable range 
for the gas mass in the core of a group (see for example, 
Fig 4 of David et al. 1995). More importantly, however, 
we must consider whether such a thin structure could sur- 
vive in the cluster environment. Following Cowie & Mc- 
Kee (1977), we have computed the evaporation rate for 
a spherical cloud embedded in the hot ICM of the Coma 
cluster. In the Coma cluster core, the classical thermal 
conduction should provide a reasonable estimate for the 
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evaporation time, which is given as 

Wap = 3.3 x l0 20 n c R 2 pc T- 5/2 (\n A/30) yr, (1) 

where n c is the embedded cloud density, R pc is its ra- 
dius in pc, Tf is the temperature of the surrounding, hot 
medium, and A is the Coulomb logarithm. For a cloud of 
radius 10 2 kpc, density 2.5 x 10~ 3 cm~ 3 , embedded in a 
hot medium of temperature 8 keV, the evaporation time 
is tevap ~ 0.1 Gyr. In 0.1 Gyr, a galaxy or group trav- 
eling at the cluster velocity dispersion could traverse 0.2 
Mpc. The evaporation time at the end of the filament 
away from the cluster center where its density is lower 
is significantly shorter. In principle, cylindrical geometry 
makes the evaporation time longer by a factor of a few. 
Also, any chaotic magnetic field could significantly reduce 
the thermal conduction efficiency and increase the evapo- 
ration time-scale. Finally, we do not necessarily need the 
evaporation time-scale to be smaller or comparable to the 
traverse time: if the amount of gas deposited into the tail 
significantly exceeded 5 x 10 11 M , it may be mostly evap- 
orated at present, with only the remnants of the group gas 
providing the observed surface brightness enhancement. 

An alternative is that the observed filamentary feature 
arises from cluster gas compressed in the potential well 
of the dark matter perturbation formed by a tidally dis- 
rupted infalling galaxy group. Merritt (1984) has shown 
that galaxy halos are tidally disrupted by the cluster grav- 
itational field. Specifically, a gravitationally bound system 
entering a cluster is tidally limited at a radius r t given by: 

r t /R c = F{r/R c )o- g /a c (2) 

where R c is the cluster core radius, a g is the group veloc- 
ity dispersion, a c is the cluster velocity dispersion, r is the 
distance from the cluster center, and F(r/R c ) is a function 
which has a minimum value of 1/2 at r — R c (see Fig 1 in 
Merritt 1984). The function F(r/R c ) decreases nearly lin- 
early with decreasing radius and hence an infalling group 
will be tidally disrupted as it approaches r = R c . The 
minimum size depends in detail on the parameters of the 
group and cluster, but for a typical group with a g — 300 
km sec -1 in the Coma cluster with a c = 1142 km sec -1 
(Biviano et al. 1996) and R c = 420 kpc (Briel et al. 1992), 
the maximum tidal radius for an infalling group is r t ~ 50 
kpc. Thus, the group will be progressively disrupted as 
it approaches the core. A trail of tidally stripped de- 
bris behind the group provides the gravitational potential 
perturbation which enhances the surface brightness. For 
isothermal gas in hydrostatic equilibrium, the gas density 
perturbation is: 

(po + Sp)/p = exp(-p J m p Aip/kT). (3) 

where po is the initial gas density, and kT is the gas 
temperature. The observed surface brightness enhance- 
ment corresponds to a gas density enhancement of about 
5%, hence the change in the gravitational potential, Aip, 
can be written as 1.05 = cxp(—pm p A(p/kT) or Aip = 
0.05kT/pm p . 

The gravitational potential on, and near, the surface of 
a long cylinder is: 

V = -G Pl \n(L/R) (4) 

where pi is the linear mass density, L is the length of the 
cylinder, and R is its radius. If we equate the potential 



on the surface of the cylinder to the change in the gravita- 
tional potential required to produce the observed surface 
brightness enhancement, then we can solve for the linear 
mass density in the cylinder. This provides an expression 
for the added mass in the cylinder: 

A(pL _ O.ObkTL 
91 ~ G\n(L/R) ~ pm p G\n{L/R) ( ' 

For the observed properties of the filament, and assuming 
the gas temperature is that of the cluster, we find a mass 
of M = 6 x 10 12 M Q . At most, the mass could be 50% 
smaller if the temperature of the gas in the filament is at 
its measured lower limit. The surface brightness enhance- 
ment also may be underestimated by a factor of a few, 
since the filament may not be located in the densest part 
of the projected gas distribution in which case the required 
density contrast is larger. Also, the 10% enhancement in 
the surface brightness corresponds to more than 10% en- 
hancement in the volume emissivity, because the filament 
is narrow and projected onto the main cluster emission. 
Thus, the measured mass in the perturbation is approxi- 
mately 6 x 10 12 Mq and is uncertain by a factor of about 
2. This mass is well within the range of that in the central 
region of a group (e.g. David et al. 1995). 

We have discussed the filament in some detail, but the 
remaining emission is also of interest. The top-right panel 
in Fig 1 compares the "raw" smoothed image of Coma 
[top-left) with the smoothed image after subtraction of 
the detected small scale structure {bottom panels). The 
removal of the two central X-ray peaks around the galaxies 
NGC 4889 and NGC 4874, as well as the filament, leaves a 
smooth and regular-looking X-ray surface brightness dis- 
tribution. A relatively small surface brightness perturba- 
tion can produce marked irregularity in the surface bright- 
ness distribution. Therefore, the small scale structure can 
significantly influence the overall appearance of the surface 
brightness distribution of the cluster. Numerical simula- 
tions must retain sufficient resolution, if they are to fully 
model shapes of the X-ray cluster images. 

4. CONCLUSIONS 

We report the detection of a long, linear filamentary fea- 
ture extending over approximately 1 Mpc from the Coma 
cluster center toward NGC 4911. We explore two possi- 
bilities for the origin of this filament. First, we show that 
the filament could arise from ram pressure stripped gas 
from a group passing through Coma, especially if mag- 
netic fields significantly reduce the conduction efficiency 
and lengthen the evaporation time. Second, we note that 
tidally stripped dark matter from a merging group could 
produce a perturbation that would give rise to the ob- 
served structure. 

Increases in sensitivity and new analysis techniques have 
continued to demonstrate the complexity of the cluster 
environment. The Coma cluster, long thought to be the 
prototype for a relaxed cluster, is in fact a remarkably 
complex system with structure on many different scales 
seen both optically and in X-rays. 

We thank L. David for helpful discussions and acknowl- 
edge support from NAS8-39073. 
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